Influence of vector interactions on the hadron-quark/gluon phase transition 
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The hadron-quark/gluon phase transition is studied in the two- phase model. As a further study of 
our previous work, both the isoscalar and isovector vector interactions are included in the Polyakov 
loop modified Nambu-Jona-Lasinio model (PNJL) for the quark phase. The relevance of the ex- 
change ( Fock ) terms is stressed and suitably accounted for. The calculation shows that the isovector 
vector interaction delays the phase transition to higher densities and the range of the mixed phase 
correspondingly shrinks. Meanwhile the asymmetry parameter of quark matter in the mixed phase 
decreases with the strengthening of this interaction channel. This leads to some possible observation 
signals being weakened, although still present. We show that these can be rather general effects of 
a repulsion in the quark phase due to the symmetry energy. This is also confirmed by a simpler cal- 
culation with the MIT-Bag model. However, the asymmetry parameter of quark matter is slightly 
enhanced with the inclusion of the isoscalar vector interaction, but the phase transition will be 
moved to higher densities. The largest uncertainty on the phase transition lies in the undetermined 
coupling constants of the vector interactions. In this respect new data on the mixed phase obtained 
from Heavy Ion Collisions at Intermediate Energies appear very important. 

PACS numbers: 12.38.Mh, 25.75.Nq 



I. INTRODUCTION 

The phase transition from nuclear matter to quark- 
gluon matter is one of the most concerned topics in mod- 
ern physics related to heavy-ion collision experiments and 
compact stars. As a principle tool, lattice QCD provides 
us a framework to investigate non-pcrturbative phenom- 
ena, such as confinement and quark-gluon plasma forma- 
tion at finite temperature and vanishing (small) chemical 
potential [1-7]. However, lattice calculations suffer the 
sign problem at large chemical potential. To evade this 
problem several approximation methods have been pro- 
posed [8-12], but the validity of the result at fiq/T > 1 
still should be taken with care [13]. On the other hand, 
to give a complete description of QCD phase diagram 
some phcnomcnological effective models [14-24] have 
been also developed. Among these models, the Nambu- 
Jona-Lasinio (NJL) type models [19-24], especially those 
coupled with Polyakov loop (PNJL) [28-35] are predom- 
inant, offering a simple illustration of chiral symmetry 
breaking and restoration, as well as (de) confinement ef- 
fect. 

The lattice QCD and (P)NJL type models are based 
on the degrees of freedom of quarks and giuons. Re- 
cently, the two-phase model with both hadron and quark 
degrees of freedom, widely used in the description of the 
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phase transition in neutron star matter under the weak 
equilibrium (e.g., [36-45] ), has also been taken to study 
the phase transition related to heavy-ion collisions [46- 
54] , particularly the phase transition in asymmetric mat- 
ter. The latter is possible to be probed in the planned 
facilities, such as FAIR at GSI-Darmstadt and NICA at 
JINR-Dubna [47, 49-54]. In these studies, only the scalar 
interacting channel was considered for quark matter, and 
this channel interaction is responsible for the dynamical 
masses of quarks. The isoscalar- vector channel interac- 
tion was also included to study the properties of quark 
matter [30, 55] or the hadron-quark phase transition in 
neutron star matter [56, 57]. 

However, up to date, the vector interacting channels 
(including both the isoscalar vector and isovector vec- 
tor interaction) have not been considered in describing 
quark matter, in the context of heavy-ion collisions, in 
the two-phase model. The inclusion of these vector in- 
teractions will modify the quark pressure and chemical 
potentials. In particular, the isovector-vector interact- 
ing channel contributes to u, d quark flavors differently. 
Correspondingly, the onset densities of the quark phase 
are possibly modified and some observation signals of 
hadron-quark phase transition in asymmetric matter may 
be infiuenced. Therefore, it is important to study the 
effect of the vector channel interactions on the hadron- 
quark phase transition, and related information deduced 
from new data in heavy ion collisions at intermediate en- 
ergies. 

The paper is organized as follows. In Section II, we 
describe briefly the modified effective quark model and 
give the relevant formulas with the newly included isovec- 
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tor and isoscalar vector interactions. In Section III, we 
present the numerical results and discuss the influence of 
vector channel interactions on the hadron-quark phase 
transition of dense asymmetric matter. For the isovector 
part we include also some results with an Isospin-MIT- 
Bag model, just to stress the physics behind the sym- 
metry terms in the quark sector, not dependent on the 
models. Finally, a summary is given in Section IV. 



II. 



THE MODEL 



In the two-phase model hadron matter is described by 
the nonlinear Walecka model, and quark matter is de- 
scribed by the PNJL model with the newly added vector 
interactions. In the mixed phase, the pure hadronic phase 
and quark phase are connected to each other through 
the Gibbs conditions with the thermal, chemical and me- 
chanical equilibria, based on baryon number and isospin 
conservation in the strong interacting process. 

The Relativistic Mean Field (RMF) approach will be 
taken to describe the properties of hadronic matter. 
This model can provide an excellent description of nu- 
clear matter and finite nuclei. The exchanged mesons 
in this model include the isoscalar-scalar meson (a), 
isoscalar- vector meson (w), isovector- vector meson (p) 
and isovector-scalar meson (6). This is called the Non 
Linear-p, S (NLpS) effective interaction. For details, see 
Refs. [41, 47-50] and references therein. 

For the quark phase, we adopt an extension of the 
two-flavor NJL model to include the Polyakov loop 
contribution[28, 29]. The Lagrangian is given by 



q{i-f''D^-mo)q + Ga 



{qqf + [qil^fqf 



-Gi. 



{qfqf + {qi-fsq)" 



[q-f^'Tqf + {q-t-.YTqf 



{qi^'qf + [qiil'^qf 



-UmA]MAlT) (1) 

where q denotes the quark fields with two flavors, u and 
d, and three colors; toq = diag(m„, m^) in flavor space. 
The covariant derivative is defined by 



On - iA„ 



with the background gluon field A^^ = (5^,o^o supposed 
constant and uniform. The temperature-dependent 
Polyakov effective potential, $[A], T), is a func- 

tion of the Polyakov loop $[^4] and its hermitian conju- 
gate $[A]. In some analogy with the nonlinear Walecka 
model, in this study we also try to include in the NJL 
term the isoscalar -vector and isovector-vector interac- 
tion channels as given in Ref. [57]. 

Here some considerations are in order. In our pre- 
vious work, Refs. [53, 54], where only the scalar inter- 
action channels were considered, calculations were per- 
formed within a relativistic mean field approximation 
which essentially corresponds to the Hartree approxima- 
tion. However, the inclusion of the Fock (exchange) terms 



would be desirable. In fact, these terms originate from 
the correlations due to the Fermi-Dirac statistics, there- 
fore they are related to a genuine quantum effect, which 
in general cannot be neglected when studying a many- 
body system. Including the Fock terms, the whole variety 
of processes in quark dynamics arising from the fermionic 
intrinsic degrees of freedom ( spin, flavor and color ), are 
automatically accounted for. Indeed, even starting from 
an effective Lagrangian containing only scalar channels, 
the exchange terms naturally yield contributions in the 
vector channels, as we will show in the following. 

A quantity of interest in the study of quark dynamics is 
the statistical average of the canonical energy-momentum 
density tensor, <: T^i/(x) :>, from which thermodynam- 
ical quantities, such as the pressure, can be derived. For 
the considered Lagrangian, the interaction part of this 
quantity, <: T^j^{x) :>, depends on the statistical average 
of the product of four quark fields, that in the Hartree- 
Fock ( HF ) approximation, can be written as: 

<: qa{x)qa{x)q^{x)q^{x) :> 

=<: qa{x)qa{x) :><: q^{x)q^{x) :> 

- <: qa{x)q^{x) :><: q^{x)qa{x) :>, (2) 

where the brackets denote statistical averaging and the 
colons denote normal ordering. It is useful to define the 
matrix [F(x)]ai3: 



<: qfi{x)qa{x) :> 



where a and /3 are triple indices for spin, isospin (fiavor) 
and color. 

As far as one is concerned with the equilibrium prop- 
erties of isotropic and non-colored quark matter it is suf- 
ficient to consider only the scalar and vector channels 
( isoscalar and isovector ). Then the matrix F{x,p) can 
be decomposed as: 



F{x) = F{x)+^^F''{x) 

+f-B{x)+-f,,T-Bf'{x) 



(3) 



It should be noticed that this matrix is related to the 
various densities characterizing the system. Indeed the 
scalar and current isoscalar densities are given by: 

ps{x) = <: q{x)q{x) :>= TrF{x) ^ 8 F{x) , 



j''{x) =<: q{x)-f''q{x) :>:= Tr-ff^Fix) ^SNcFf'ix) , 
while the isovector counterparts are given by: 

P3{x)^<: q{x)T3q{x) :>= TrTsFix) =8N,B3{x), 

j^ix) = <: qix)Tsj'^qix) :>= Tnf^T^Fix) =8N,B^ix) , 

where the traces are taken over spin, fiavor and color 
indices. Then the interaction part of the energy- 
momentum density tensor, in the HF approximation, 
reads: 
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<:TW(a;) :> = Tli\x)Hartre 



2G„Tr{F{x)F{x) + {i-i5)TF{x) ■ {ij5)rF{x)) 



-2G^Tr{jxF{x)j^F{x)+j5lxF{xh5l^F{x)) + 2GsTr{fF{x) ■ tF{x) + {i^^)F{x){i-i^)F{x)) 



^2GpTr{fjxF{x) ■ f-f^F{x) + fj5lxF{x) ■ Tj5l^F{x)) 



(4) 



The exchange terms that appear in the energy- 
momentum tensor can be evaluated exploiting the de- 
composition given in Eq. (3). After some algebra one 
realizes that the effect of the Fock terms is equivalent to 
redefine the coupling constants, as written below: 

Gcr ~ Ga + {Gcr — Gs)/12 
Gs = Gs — (Gcr — Gs)/12 
G^ = G^{1 + 1/6) + (G, + Gs)/6 + Gp/2 



Go 



Gp(l-l/6) + G„/6 . 



(5) 



Then, with the effective coupling constants given above, 
calculations can be performed as in the Hartree approx- 
imation. 

It should be remarked that all the relevant interaction 
channels in general can occur in the HF approximation, 
even if some channel is absent in the original Lagrangian. 
For instance, if the vector channels are not present, i.e., 
we take G^j = and Gp = as in our previous work, con- 
tributions to the uj channel are naturally arising from the 
exchange terms associated with the scalar channels. On 
the other hand, the vector isovector p channel gets con- 
tributions (both direct and exchange) only from vector 
channels. 

Since the critical end-point of the first order chiral 
transition appreciably depends on the strength of the 
vector channel interaction [58-60], the exchange contri- 
bution (Gct -I- Gs)/Q to the effective value of G^j could 
represent a reference value for the vector channel inter- 
action [61]. 

In the following we will adopt the choice Gs — and 
the notation G = G^ = (G^ + Gs)- The coupling con- 
stants of the vector interactions, G^ and Gp, will be taken 
as parameters, and different values will be used to inves- 
tigate their influence on the phase transition. For conve- 
nience we define r^j = Gu>/G, Tp = Gp/G. 

For the temperature dependent effective potential 
U{^, we use the parametrization given in Ref. [64] 



iY($,$,r) 

ji4 



_a(T) 
2 

+4(1.3 



-I- 5(r)ln[l - 6$$ 
+ $^) -3($$)2], (6) 



where 



a[T) = flo+fli 



(7) 



The parameters Oi, bi are fitted to the lattice QCD re- 
sults in pure gauge theory at finite temperature. In the 
equation above Tq represents the temperature where the 
Polyakov potential gives a deconfinement phase transi- 
tion in a pure gauge theory. The original value of Tq 
fitted to pure gauge lattice QCD data is 270MeV [65]. 
When fermion fields are included, the temperature Tq is 
usually rescaled to obtain a consistent result with the full 
lattice data, which give the value T'^ = 173 ± 8MeV for 
deconfiniment transition temperature [1. 2, 4]. In this 
paper the value 210MeV for Tq is adopted. 

The PNJL model is not renormalizeable, so a cut-off 
A is introduced to get finite results for three-momentum 
space integrations. For the model parameters we take the 
values A = QblMeV, G = 5.04GeV-^ m„,d = 5.5MeV, 
determined by fitting the chiral condensate, /^r and 
to their experimental values [29] . The coefficients in the 
Polyakov effective potential are listed in Table I. 



TABLE I: Parameters in Polyakov effective potential given 
in [64] 
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3.51 


-2.47 
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-1.75 



The thermodynamical-potential density of quark mat- 
ter in the mean field approximation reads 



n = W($, $, T) + G(0„ + ^df - GUPu + Pdf - Gp{pu - Pdf 



(2^ 



3(£;„ + Ed) 



(2^ 



ln(l + 3<i>e-(^'-^*)/^ + 3$e-2(^--''*)/r + e-3(£.-M*)/T) 
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ln(l + 3$e-(^'+^*)/^ + 3<i>e-2(^^>+'^*)/T ^ g-3(B.+M:)/T^ 



(8) 



where pi is the number density of quarks of flavor i, 
Ei = 2 + Mf and /i* are the corresponding energy- 
momentum dispersion relation and effective chemical po- 
tential with 

Hl = Pu-^Gu,{pu + Pd) -2Gp{pu- Pd) (9) 

p*i = pd- '2G^{pu + Pd) + ^Gp{pu - Pd) (10) 

lynamical quark masses and qu; 
coupled with the following equations 



M = mo - 2G(0„ + 0d) 
d^k M 



(27r)3 E 



(11) 



(l-n,(fc)-fi,(fc)), (12) 



where ^^(fc) and fii(fc) are modified Fermion distribution 
The dynamical quark masses and quark condensates are functions of quark and antiquark (similar to those given 

in Ref. [66] without the vector contributions) 



ni{k) 



(13) 



ni{k) 



1 + 3l>e-(^'+''?)/^ ■ 



I 



(14) 



The values of 0„ , 0^ , $ and $ are determined by mini- 
mizing the thermodynamical potential 

_ dn on on _ 

Wu'Wd^d^^M^^- 

All the thermodynamic quantities relevant to the bulk 
properties of quark matter can be obtained from 57. Par- 
ticularly, we note that the pressure and energy density 
should be zero in the vacuum. 

The number density of quarks of flavor i 



2x3 



(2^ 



(ni(fc) - ni{k)) 



(16) 



can be derived by means of the relation pi = —dfl/dpi. 
The baryon and isospin densities and the corresponding 
chemical potentials in quark phase are defined by 



Pb = ^(Pu + Pd), 



PS 



Pd, 



:{p.u + p-d), Pa = 7;iP'u - P-d) 



(17) 



(18) 



while the asymmetry parameter of quark matter is de- 
fined by 



Pb 



Pu + Pd 



(19) 



Analogous definitions hold for hadronic matter: 

Pn, P3^Pp-pn, (20) 



Pb = Pp 



H 1/ 

Pb = -j^KPp 



Pr, 



pi 



1 



{Pp Pn) 7 



Pn - Pp 
Pn + Pp 



(21) 



(22) 



where pp and p„ are the proton and the neutron densities, 
respectively. When a mixed phase of quarks and hadrons 
is considered, the Gibbs' conditions (thermal, chemical 
and mechanical equilibrium) 



Pb{Pb^PziT) = p%{pb,p^,T) 
PaiPB^Ps^T) = ^^(pb,P3,T) 
P"{pB,P3,T) = pQ{pB,p,,T), 



(23) 



should be fulfilled (a general discussion of phase tran- 
sitions in multicomponent systems can be found in 
Ref. [36]). In Eq. (23), = (1 - x)Pb + XpI is the 
total baryon density and P3 = (1 — x)pi -I- XP? is the 
total isospin density, where x is the quark fraction. In 
heavy-ion collisions, for a given initial charge asymmetry 
the global asymmetry parameter a of the mixed phase 



_ P3 _ (1 - X)pi + XP? 



Pb (1 - x)pI + XP3 ■ 



(24) 
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should be constant according to the charge conservation, 
but the asymmetry parameters a^,a'^ in the separate 
phases can vary with x- For details, one can refer to 
Refs. [49, 53, 54]. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

In this section, we display the numerical results and 
discuss the influence of isovector and isoscalar vector 
interactions on the phase diagram of the hadron-quark 
phase transition. 

In actual calculations a value of 0.2 is chosen for the 
asymmetry parameter a. We notice that in heavy-ion 
collision experiments the largest value of a, a = 0.227, is 
possibly reached in ^ssy _|_238 u collisions. 



A. The role of the isovector vector interaction 

Firstly, we only focus on the influence of the isovector 
vector interaction on the phase transition in the two- 
phase model, then we set — 0, and perform calcu- 
lations for various values of Vp, Tp ~ 0.0, 0.25, 0.5, 1.0. 
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FIG. 1: (Color online) Diagram of the hadron-quark phase 
transition in the T — pg plane for symmetric (a = 0) and 
asymmetric matter (a = 0.2) with Vp = 0.0, 0.25, 0.5, and 
1.0. The baryon density pg is expressed in units of the density 
of ordinary nuclear matter po- The lines in the left side cor- 
responding to X = represent the onset of the mixed phase, 
and those in the right side corresponding to x = 1 denote 
the beginning of the pure quark phase. The dot indicates the 
critical end point. 



We display the phase diagram in T— plane in Fig. 1. 
In agreement with the results of Refs. [53, 54], this figure 



shows that the onset density of the mixed phase in asym- 
metric matter is smaller than that of symmetric matter. 
This means that in heavy-ion collision experiments the 
phase transition could be relatively easier to be reached 
for asymmetric matter with respect to the symmetric 
case. However, the onset density appreciably depends 
on the strength of the interaction in the isovcctor-vector 
channel. It is shifted to higher densities with increas- 
ing the coupling parameter Gp. Meanwhile the density 
range of the mixed phase shrinks. Indeed, similarly to 
what is observed in the hadron sector, the isovector vec- 
tor channel yields a positive (repulsive) contribution to 
the pressure of the quark phase, connected to the corre- 
sponding symmetry energy. 
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FIG. 2: (Color online) The same of Fig. 1 but for the T- 
plane. 



Similar features are observed for the phase diagram in 
the T — plane as shown in Fig. 2. However, we notice 
that for symmetric matter, only one phase-transition line 
exists. This is due to the fact that in the symmetric case 
the curve T — hb does not depend on the quark fraction 
X, see also Refs. [53, 54]. Both figures 1 and 2 show 
that the isovector-vector channel of the quark interaction 
plays a more important role on the phase transition with 
increasing density and lowering temperature. In fact, 
at higher temperatures and/or lower densities the role 
of interactions in general becomes weaker with respect 
to the kinetic contributions. Finally, we still observe the 
occurrence of a critical end point (CEP) of the first order 
phase transition, analogously to what happens when only 
scalar channels are considered [53, 54]. 

In Fig. 3 the asymmetry parameters of hadronic and 
quark matter in the mixed phase are displayed as a 
function of the quark fraction x^ for the temperature 
T = 100 MeV . We can observe a clear Isospin Distil- 
lation effect [47, 49], i.e. the asymmetry of quark matter 
is much larger than 0.2 at the beginning of the phase 
transition and decreases with increasing the quark frac- 
tion. Whereas the asymmetry of the hadronic matter 
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keeps below 0.2 and is a slowly decreasing function of x- 
These features of the local asymmetry may lead to some 
observable effects in the hadronization during the expan- 
sion phase of heavy ion collisions, such as an inversion 
in the trend of emission of neutron rich clusters, an en- 
hancement of 7r^/7r+, /\"//\+ yield ratios in high-density 
regions, as well as an enhancement of the production of 
isospin-rich resonances and subsequent decays, for more 
details see Refs. [49, 53, 54]. These signals are possible to 
be probed in the newly planned facilities, such as FAIR 
at GSI-Darmstadt and NICA at JINR-Dubna. 
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FIG. 3: (Color online) Asymmetry of hadronic and quark 
matter in the mixed phase as a function of the quark con- 
centration for various values of the coupling constant in the 
isovector-vector interaction channel. 



pulsive term in the quark sector. In order to confirm that 
this is a general result, not depending on the different 
quark models, we present here few similar calculations 
of the hadron-quark transition, in a two-EoS approach, 
using the MIT-Bag model for the quark matter [47-50]. 

Now the results are very sensitive to the choice of the 
Bag-constant B, in particular at low temperatures and 
high baryon densities. In fact for low B- values we can 
even get a disappearing of the transition since the hadron 
pressure cannot match anymore the quark pressure, as 
discussed in detail in Ref. [50]. 

Here we choose a Bag-constant B = {l&QMeVY , 
which gives for symmetric matter a high-density mixed 
phase structure very close to the one of Fig. 1, obtained 
with the PNJL model. Of course the same NLpS Rel- 
ativistic Mean Field interaction is used for the hadron 
sector. 

A quark isovector-vector (p-like) term is introduced by 
a naive application of a constituent quark model of the 
nucleons, i.e., just reducing of a factor 3 the Nucleon- 
p coupling constant of the hadron part. As a conse- 
quence we get straightforward corrections to the quark 
pressure and chemical potentials with respect to the sim- 
ple relativistic Fermi gas values of the MIT-Bag model 
ipQ{F),fi^iF),fi,iF)y. 



M« = ^^u{F) + '-^{pu ~ Pd), 



ml 

f^d = f^diF) - ^{Pu - Pd), 
mi 



(25) 



However we also see that the strengthening of the 
coupling constant Gp reduces the distillation effect and 
consequently the asymmetry parameter a'^ , which may 
weaken the observational signals of the phase transition, 
although always present. The uncertainty is that the 
relevant p-coupling constant cannot be unambiguously 
determined from lattice-QCD calculations. So the ob- 
servation or not of the related signals in experiments can 
provide some hints on this aspect. 

In any case we note that both isospin effects, earlier 
density transition and isospin distillation are still there 
even with large values of Vp. They appear as general, 
quite robust effects of asymmetric matter, not much af- 
fected by the introduction of a symmetry interaction 
term in the quark sector. In order to confirm all that, in 
the next paragraph we present other results obtained in 
a much simpler Iso-MIT-Bag model. 



B. Few results with the MIT-Bag model 

As already remarked the isospin effects on the mixed 
phase (boundaries and asymmetries in the two phases) 
are naturally related to the presence of a symmetry re- 



with Qp^q = gp^N/'i and mp is the p-meson mass. The 
coupling choice is fixed by the NLpS parametrization of 
the hadron sector, 



Ip.N 



yp,N 



3.15/m2 



see the detailed Appendix A of Ref. [50]. We note 
that such simple insertion of isovector-vector terms in 
quark phase (Iso-MIT-Bag model, in the following named 
ISOE results) can be particularly justified at high 
baryon densities and chemical potentials, where we ex- 
pect a more relevant role of the hadronic degrees of 
freedom. In any case here our aim is to show general 
quark symmetry energy effects on the hadron-quark tran- 
sition at low temperatures in isospin asymmetric matter, 
that confirm the results of the previous section with the 
Hadron-PNJL approach. 

The Fig. 4 is the T ~ pg phase diagram for symmetric 
and a = 0.2 isospin asymmetric matter, corresponding 
to the previous Fig. 1 obtained with the PNJL model in 
the quark sector, with various weights of the p coupling. 
We clearly sec that at temperatures below 40 MeV and 
at high baryon densities the curves are very similar, in 
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particular in the choice in Fig. 1 of a rp ratio equal to 0.5, 
close to the evaluation around 0.7 we use in the MIT-Bag 
model following the constituent quark picture. 

The same comment is valid for the behavior of the 
asymmetry parameters in the mixed phase, shown in the 
ISOE calculation in the Fig. 5, to compare to the previous 
Fig. 3. 
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FIG. 4: (Color online) iVLptJ-MITBag results: T-pg plane of 
the hadron-quark phase transition for symmetric (a = 0) and 
asymmetric matter (a = 0.2) with the ISOE choice for the 
p-coupling. The curves in the left side with x — represent 
the beginning of the mixed phase, and those in the right side 
with X — 1 mean the beginning of the pure quark phase. 
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FIG. 5: (Color online) A^Lp^-MITBag results: Asymmetry of 
hadronic and quark matter in the mixed phase as functions 
of quark concentration without and with the ISOE choice for 
the isovector vector interaction coupling. 



In general we can say that all the isospin influence on 
the hadron-quark transition is still present when we in- 
troduce vector-isovector terms in the quark phase, even 
with relatively large weights. Of course the interaction 
symmetry repulsion in the quark sector will reduce both 
effects, earlier transition densities and isospin distillation 
in the mixed phase, but the possibility of related obser- 
vations appears still there. 



C. The role of the isoscalar vector interaction 

We discuss now the results obtained when the 
isoscalar-vector interaction channel in the quark sector 
is turned on. We choose the value of 0.2 for the ratio 
= Gu/G. This value is close to the contribution to 
this channel from the exchange terms of the scalar chan- 
nels, see Eqs. (5). 
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FIG. 6: (Color online) The same as Fig. 1 except for the 
inclusion of isoscalar vector interaction with r^i = 0.2 



In Figs. 6 and 7 are displayed the diagrams of the 
hadron-quark phase transition in the T — ps and T — plb 
planes respectively. Compared to Figs. 1 and 2 the 
phase-transition curves are significantly moved toward 
higher values of density /chemical potential. This can be 
explained in terms of the repulsive contribution of the 
isoscalar-vector channel to the quark energy and, as a 
consequence, to the chemical potential ( see Eqs. 9 and 
10 ). More specifically, the relevant quantity in the ki- 
netic contribution to the thermodinamical potential, be- 
sides the temperature, is the effective chemical poten- 
tial, p*. This quantity is appreciably smaller than the 
full chemical potential. Then, higher values of the lat- 
ter quantity are necessary to fulfil the Gibbs conditions. 
Moreover, we obeserve that for low values of temperature 
the quark chemical potential attains values of the same 
order of magnitude as the cut-off A. However, this does 
not give rise to inconsistency, since the momentum scale 
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is set by the Fermi momentum, which is determined by 
the effective chemical potential. Also the CEP moves to 
higher values of the density /chemical potential although 
slightly, whereas the corresponding critical temperature 
almost keeps the same value. Moreover, in asymmetric 
matter the region of coexistence of the two phases is more 
extended when the isoscalar-vcctor interaction channel is 
included in the quark sector. 
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inclusion of the repulsive vector interaction in the quark 
sector affects significantly the possibility to observe the 
first order hadron-quark phase transition, especially in 
the low temperature-finite chemical potential region. In- 
deed, with = 0.2, the onset transition density, even in 
asymmetric matter, is already shifted to density values 
larger than 5.5 po- Thus the study of this transition and 
of the associated phase diagram may give information 
on the properties of the interaction between quarks and 
corresponding coupling constants. From the definition 
of the effective coupling constants, Eqs.(5), one can see 
that a considerable positive contribution to G^n equal to 
0.17 G, already comes from the exchange terms of the 
scalar channels. On the other hand, the observation of 
such a transition in relativistic heavy ion collisions at 
intermediate energies, where moderate density and tem- 
perature values are reached, would point to smaller value 
of Vi^. Thus such an evidence could be taken as an indi- 
cation of the presence of the vector channels already in 
the original quark Lagrangian, with a negative Guj and a 
positive Gp. We note that the sign of Gp is dictated by 
physical reasons: it must be positive in order to obtain a 
positive symmetry energy. 



FIG. 7: (Color online) The same as Fig. 2 except for the 
inclusion of isoscalar vector interaction with r^i = 0.2 

Finally, comparing Fig. 8 with Fig. 3 one can observe 
that the inclusion of the isoscalar-vector channel leads 
to a small increase of the asymmetry parameter of the 
quark matter in the mixed phase. The Isospin Distilla- 
tion effect may be strengthened by this interaction chan- 
nel, but the onset of the mixed phase is shifted toward 
higher densities. 
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FIG. 8: (Color online) The same as Fig. 3 except for the 
inclusion of isoscalar vector interaction with r^i = 0.2 



From the results discussed above it is clear that the 



IV. SUMMARY 

We have studied the hadron-quark/gluon phase transi- 
tion in the two-phase model with the newly added isovec- 
tor and isoscalar vector interactions for quark matter. 
We stress the presence of vector terms in the quark Equa- 
tion of State just due to exchange contributions from 
scalar fields in the effective PNJL Lagrangian. 

The consideration of the isovector vector interaction 
splits the effective quark chemical potential of u and d 
quarks. The calculations show that the phase transition 
densities are delayed to higher values and the ranges of 
the mixed phases shrink with the increase of the coupling 
constant Gp. Meanwhile the asymmetry parameter a'^ 
at small x is reduced for a larger Gp. 

Furthermore, with the inclusion of the isoscalar vector 
interaction, the whole phase diagrams move to higher 
baryon densities/chemical potentials, but the tempera- 
tures of the CEP almost keep unchanged. 

In our previous study, we have proposed some possi- 
bly observable signals of the hadron-quark/gluon phase 
transition for asymmetric matter in heavy-ion collision 
experiments. These signals are possibly weakened when 
the isovector vector interaction is included, whereas they 
are slightly strengthened by the isoscalar vector interac- 
tion. However, the main problem induced by the isoscalar 
vector interaction is that the onset densities of the mixed 
phase are moved to higher densities. Most uncertainty 
lies in the relevant vector couplings. The planned ex- 
periments with FAIR at GSI-Darmstadt and NICA at 
JINR-Dubna are expected to provide some hints on the 
related study. 
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